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PREFACE 

The work reported herein was performed for the National Aeronautics and 
Space Administrations's Marshall Space Flight Center under contract 
NAS8-33191, Supplemental Agreements 2 and 3. This effort involved the 
development of a computer program to perform Statistical Energy Analysis. 
This volume constitutes the final deliv2rable item under the contract. A 
card deck of the computer program was forwarded previously by MDAC letter 
A3-130-GWJ-1283. 

Programming was accomplished by S. J. Nygaard, McDonnell Douglas Automation 
Company - Huntington Beach. 
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Sect i or 1 
INTRODUCTION 

Significant hiqh frequency random vibration environments are qenerated 
durinq the operation of aerospace vehicles within the atmosphere. To 
achieve optimum vehicle desiqn, vibration and acoustic criteria are 
developed early in the vehicle development proqram and are updated 
periodically as the desiqn matures. 

Hiqh frequency random vibration response prediction does not lend itself 
well to classical structural dynamics and such predictions are usually 
made by extrapolation from existinq data banks. This method qives 
excellent results when similar structures are involved. However, as 
similarity decreases and associated extrapolations become larqe, uncer- 
tainty over accuracy also increases. 

The efforts of numerous investiqators (Maidanik, Lyons, et al.) have 
examined a more qeneral hiqh frequency random vibration analysis 
approach - the so-called Statistical Enerqy Analysis method. The SEA 
method is able to accomplish hiqh frequency prediction of arbitrary 
structural configurations and is therefore a significant improvement over 
extrapolation methods when little or no previous data exist. The SEA 
method also represents a qreat Improvement over normal mode methods for 
hiqh frequency random vibration prediction because of the qreatly reduced 
computational complexity associated with the more qeneral SEA model 
elements and the attendant Improvement In analysis turnaround time. 

SEA has been developed for complex structures by MDAC under contract to 
MSFC (Ref, 1 and 2 ) over the oast three years. The past year's effort 
has created a qeneral SEA computer proqram which is described In this 
manual. The manual contains (1) a summary of SEA theory, (?) example 
problems of SEA proqram application, (3) a description of the computer 
proqram, and (4) a complete proqram listing (Appendices). 
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Section 2 

A SUMMARY OF STATISTICAL ENERGY ANALYSIS PRINCIPLES 


Statistical Energy Analysis (SEA) is a powerful tool for estimating the 
high frequency vibration spectra of complex systems. The analysis method 
is based on the estimation of the power flow between idealized gross 
elements of a vibrating system. The method is statistical in that 
averaging assumptions are made with regard to distribution of energy 
within an element, distribution of resonant modes, and the coupling 
between elements. T hese assumptions greatly simplify the computational 
complexity associated with normal mode methods. These same assumptions 
impose the limitation that pcint response predictions cannot be made. 

The assumptions on which the method rests and their implications can be 
quite rigorously stated as follows: 

1. The total vibrating system can be partitioned into SEA elements 
(with suitable boundary conditions) whose modes approximate the 
modes of the original vibrating system. 

2. The modes of the elements of a system contain all of the vibratory 
energy of the system. 

3. The energy in one frequency band of a system element is equally 
distributed among the modes of that element occurring in the 
frequency band. 

4. Only modes occurring within the same frequency band are coupled. 

5. For two coupled elements, all of the modes occurring In one of the 
elements in one frequency band are equally coupled to each mode 
occurring in the same frequency band in the other element. 


Assumption 1 contains the fundamental existence basis for SEA: the con- 
cept of parti tlonablllty. This concept Implies that a coupled vibrating 
system with system modes can be approximated by two or more separately 
Idealized vibrating elements, each with Its own independent mode set. 
These sets are coupled only In the sense of having power flow to and 
from each set across the partition boundary (later referred to as the 
"joint"). The approximation to this model exists in most structures 
having reflective boundaries in the higher frequencies. For example, a 
skin/stringer structure has higher order skin panel modes that are nearly 
the same frequency and shape as an ideally supported panel because the 
stringer is a comparatively massive boundary causing reflection of 
flexural waves from the skin panel. An SEA plate element could logically 
be equal to the panel area bounded by stringers or frames. Such elements 
will then have to be coupled with joint elements in order to develop an 
SEA model which emulates the vibratory power flow of the real structure. 

Assumption 3 is the most important simplifying assumption of SEA because 
it eliminates the necessity to calculate generalized modal forces and 
responses. The conditions implicit in this assumption are usually ap- 
proximated by the higher order modes of a structure in a reasonable 
bandwidth, say 1/3 octave. One-third octave bands represent a reasonable 
compromise between the necessity to get a fairly large number of modes 
(>20) in the band for good statistics and the necessity to have some 
frequency response resolution in the vibration prediction. The number of 
modes in a unit bandwidth can be estimated for simple structural forms 
(such as beams, plates, etc.) using algebraic expressions for modal den- 
sity such as those given in Section 4 of this report. Estimation of 
modal density in this way is a considerable simplification over normal 
mode methods. 

Given SEA elements with the properties described above it is now neces- 
sary to join them to permit power flow between the modes of one element 
and the modes of another. This is done with a parameter called the 
coupling loss factor $ and leads to assumptions 4 and 5. Assumption 4 
is directly linked to assumption 2 and the further assumption of a linear 
process. Assumption 5 follows directly from assumption 3 as part of the 
simplification associated with a statistical rather than explicit de- 
scription of modes. 
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With these properties and assumptions we can now write the SEA equation 
systems as follows: 
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where 

M is the number of SEA elements 


-N . <b. • 
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Ni = number of modes resonant in element i 
ni = element i loss factor 

4>i j » power transfer coefficient for coupling between modes in 
elements i and j 

u> » center frequency of bandwidth 

<a7S 

^i * m i — 5 — ■ total energy of element 1, <a a > being the mean square 
“ acceleration 

Si = external acoustic or mechanical excitation in the bandwidth of 
interest 


Note that the matrix a is square but not symmetric. The lack of symmetry 
arises from the nature of the term -N^j. Using the first two rows as 
an example the power balance equations are 
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Note that the "symmetric" positions and actually carry 

the number of modes belonging to the row number only. This unsymmetric 

form preserves the power flow both to and from each element. 

Each equation in the matrix states the following relationship for each 
element. The net power flow into element i (S^) equals the difference 
between (1) the power flow dissipated within the element (un^E^) plus 
the power flow lost to other elements across the joints 

M 

E i l Vlk 

k-1 

and (?) the power flow added to element i from all other elements 
M 

N i I E kV 1 * k 

k«l 


It must be remembered that the coupling term <t>. k has nothing to do 
with coupling modes; it only relates the fractional amount of energy 
resident in the modes of element i that flows to the modes of element k. 

For all but the simplest of systems even the SEA equations can be labori- 
ous to evaluate by hand. The computer program described in the following 
sections performs all of the computations necessary to evaluate terms and 
solves the above system of equations for element energies and prints the 
results as vibration PSD or RMS levels for all elements. This program 
relieves the analyst of the necessity to have an extensive knowledge of 
how to compute SEA parameters and permits the use of SEA model sizes that 
would otherwise be intractable. The program is described in detail in 
Section ?; some SEA modeling examples and corresponding program inputs 
are described in Section 3. 



Section 3 

EXAMPLES OF SEA PROGRAM APPLICATIONS 


An SEA computer program has been developed that performs the organization 
of specific problem solutions for up to a 20-element system. A number of 
different elanent and joint types are available to describe a wide vari- 
ety of structural forms in terms r ~ an SEA model. Random acoustic or 
mechanical excitation can be applied in 1/3 octave bands to any arbitrary 
number of elements within a given model. The resulting equations are 
then solved giving the vibration response spectrum far each element. 

This section illustrates the use of SEA and the computer program with 
specific examples. The basic analysis procedure Is divided into three 
steps: 

1. Idealization 

2 . Parameter Generation 

3. Problem Solution 

The idealization step must be performed entirely by the user as it 
consists of modeling the physical structure In terms of available SEA 
program elements - a conceptual process. The second step Is one of 
simply providing the proper data to the program which then carries out 
the third step. These steps will be Illustrated below with specific 
examples, and the SEA program elements will be discussed in detail. 

The idealization step is by far the most crucial step in the process. It 
is here that the art of engineering judgement must reach a well developed 
state, balancing the realities of the structural article to be analyzed 
with the capabilities and assumptions implicit in the SEA process to 
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obtain a useful engineering solution. Consider the model shown In Figure 
3-1. This Is the basic form of all SEA models. It consists of elements, 
denoted as boxes, which may be plates, beams, etc., and joints denoted by 
connecting lines which correspond to the physical Interfaces at the se- 
lected partitions. 

The model shows the articles Included In each SEA element and the con- 
nection relationships between or among the elements. To gain a better 
understanding of what these elements and connections are It Is helpful to 
now explain In detail the various types available In the computer program. 
Each of the elements Is made up of one or more sub-elements, the first of 
which Is the main sub-element. The sub-element system provides a conve- 
nient way to compute and Include the modal density and mass properties of 
Individual structural pieces that make up the element on a piece-by-piece 
basis. However, the sub-element with the most Important property being 
modelled In that SEA element should always be the first or main sub- 
element. The joint properties must also be consistent with this 
sub-element as it is the only one that can be coupled to other SEA main 
elements. An example is given by element 1 of the model. This element 
has structure elements exposed to an acoustic field and those that are 
not. The most Important property is the reception of acoustic excitation 
and transmission of vibration through its boundaries to other main 
elements. That portion must therefore be sub-element 1. 

It is often the case that structures are made of different materials; 
this condition Is taken into account when the program sets up the so- 
lution. The program automatically recomputes properties (thickness, 
density, etc.) to match the elastic modulus of element 1, sub-element 1. 
Furthermore, only sub-element 1 of each element can receive acoustic 
excitation; the others contribute only to the modal density and to the 
element mass. Sub-element 1 will always be used for the structure 
element that Is the main piece of the given SEA element because all other 
properties except mass and modal density will be those of this sub- 
element. 

The types of sub-elements available are beams, plates, cylinders, mem- 
brane, and reverberant room. In general, beams, plates and cylinders can 
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Figure 3-1. SEA Model Elements for Acoustic Test Configuration 



be freely mixed in describing a structure. The membrane sub-element 
should only be used as an ancillary sub-element (i.e., sub-element number 
> 1). The room sub-elemer.t requires special treatment in that the loss 
factor must be developed from the reverberation time (n » 2.2/fTgQ) and 
the answers must be converted from q 2 to pressure 2 usinq enerqy 
density relationship. The room sub-element must always be the only 
sub-element in that particular element. 

Given these SEA elements, joint properties must now be developed to 
describe the oower flow from the modes of one element to another. 

The proqram qives the user a choice of four types: (1) plate to plate, 

( 2 ) beam to plate, (3) f jlted joint, and (4) plate to acoustic. Each 
joint has two ends for accounting purposes and the word order used de- 
scribes the A and B ends respectively. For example, in joint (2) the 
beam is always at the A end and plate is always at the B end. The joint 
loss factors are from the literature (cited in References 1 and 2) except 
fo" the bolted joint. This joint is a plate-to-plate joint with an addi- 
tional insertion loss to account for internal losses in the joint due to 
fastener effects. The insertion loss is a load sheet input which 
provides to the user a more qeneral purpose alternative to joints (1) and 
( 2 ). Values of the inertion loss parameter for various fastener 
arrangements are not well defined, however, and must be the topic of 
continuing research. In the case of the Materials Experiment Assembly 
(MEA) analyzed durinq the last phase of this study, the empirically 
determined insertion loss factor was approximately 10 for each bolted 
joint. 

With the basics of the program now given, two examples will be shown as a 
quide for proqram use. The first example will be a seqment of skin 
strinqer structure exposed to an acoustic field with an equipment panel 
on the opposite side. The second example will be a simplified MEA 
analysis. 

The first case can be Idealized with two SEA elements as shown in Fiqure 
3-2. The relevant parameters are shown in the followinq table: 
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.040" 

18" x 72" 

10" O.C. 

1-1/2" high x 1" wide with 3/4 
flanges,. 063 thick (full hat section) 
400 in 2 of .063 thick aluminum 
1" honeycomb with .020 face sheets 
17 x 30 with 16 lb of small equip- 
ment items mounted on its surface. 
Panel riveted to frames at four 
places (U-shaped channels) .063" 
thick x 1-1/2" high 

The SEA element 1 shown in Figure 3-2 will consist of all structural 
elements in the table except the panel and its mounted equipment which is 
SEA element 2. From this information an elementary SEA analysis can be 
made. The load sheet entries are determined as follows and are shown in 
Figures 3-3 through 3-10 in the proper sequence. 


Skin Thickness 
Segment Dimensions 
Stringer Spacing 
Stringer Dimensions 
Internal Frames 
Panel 




Element 1 

Element 2 

Acoustics. 

1 

i External Skin j 
J _and Stri ngers , 1 
Internal Frames 


Equipment Panel, 

^ixcitati^^ ' 


Mounted Components 


Figure 3-2. Two-Element SEA Model Example 


Figure 3-3 shows the header card which basically describes the problem as 
consisting of two elements and eleven 1/3 octave bands from 250 to 2500 
Hz and that the output will be a vibration PSD. Figure 3-4 describes 
element 1 as consisting of two sub-elements, the first of which is ex- 
posed to acoustic excitation, and describes the damping vs. frequency 
curve. Figure 3-5 shows the 1/3 octave sound pressure levels; Figures 
3-6 and 3-7 describe the sub-element properties in detail. 
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Figure 3-3. Card 1, Example Case 



LOADSHEET (2): SEA ELEMENT PROPERTIES CARD 

o One required for each element (2 s E < 20) 
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LOADSHEET (3): EXCITATION SPECTRUM CARD 



t 


Figure 3-5. Card 3, Example Case 



LOADSHEET (4): SUB-ELEMENT PROPERTIES, CARDS 1 AND 2 
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LOADSHEET (4): SUB-ELEMENT PROPERTIES, CARDS 1 AND 2 

° Must be consistent with Loadsheet 2, CC 3-4 (SE i 1) 
° Two cards (records) per sub-element 
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Sub-element 1 consists of the skin and stringer elements. The skin is 
used as the principal property because It Is the major element being 
excited by the acoustic field. The stringers are added as a smeared mass 
because their Internal resonant frequencies are very high and therefore 
would be expected only to load the skin in the frequency range of in- 
terest. Sub-element 2 (Figure 3-7) represents the properties of the 
internal frames and channels which mount the panel. 

Figures 3-b and 3-9 show the loadsheets for the inside equipment panel. 
Figure 3-8 defines the panel as Element 2, with one sub-element and con- 
stant damping vs. frequency. Figure 3-9 has the panel properties which 
have been equivalenced to an isotropic plate because the panel is a 
composite structure. Note that both the thickness and the density have 
been changed to be consistent. Any other equivalence could also be used, 
e.g. leave the thickness = 1,0 and the density and Young's modulus will 
be recalculated using 


1 r T 

T L e iTi * Peff 
i 

12EO-Y 2 ) lit 

and ^eff * 5 — - — 

eTT bt 3 

There is no requirement that sub-elements be of similar materials, but 
all plate, beam and cylinder sub-element properties must always be iso- 
tropic equivalents. 

The weight of the components is treated as a non-structural mass. If 
half or more of the panel area were covered with these components, the F 
should be changed to a T as shown on the sheet to account for the reduc- 
tion of modal density associated with stiffening of the panel by the 
components as reported in Reference 2. 

Figure 3-10 describes the joint properties between the skin and the 
plate. Note that an added Insertion loss factor of 2 Is used to account 
for rivet effects. 
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LOADSHEET (4): SUB-ELEMENT PROPERTIES, CARDS 1 AND 2 

• Must be consistent with Loadsheet 2, CC 3-4 (SE z 1) 
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Figure 3-9 . Card 7, Example Case 



LOADSHEET (5): JOINT PROPERTIES 
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59*68 Real Energy reduction factor (BJ only) 



•? 


These data are then assembled as a file which the program reads as an 
input. 

The corresponding program output is shown in Figure 3-11 and consists of 
a labeled list of the Input data, a display of the element modal densi- 
ties, and a table of vibration responses. The adequacy of the SEA model 
with regard to assumptions 4 and 5 can be checked using the modal density 
tables of the figure. For example, element 1 contains more than 20 modes 
per 1/3 octave over the entire analysis range, whereas element 2 has far 
fewer. Since frequency response will be smoother with more modes per 
band. It is expected that the panel prediction will be poorer than the 
skin predictions. Specifically, the actual panel response may have some 
peaks which exceed the SEA prediction. 

The SEA response prediction is shown in Table 3-1. The levels for ele- 
ment 1 are high compared to the criteria published in Reference 3, but 
one must remember that this estimate includes a space average over the 
skin. The stringer and frame vibration levels which are Inputs to the 
panel are a factor of ~1C0 less. 

Simple changes can be made to improve the prediction. For example, an 
internal acoustic field can be put on the panel by also placing an A In 
card column 5 of Figure 3-8 and entering the appropriate table immedi- 
ately after, as shown in Figure 3-12. The results of this modification 
are shown in Table 3-2. As can be seen, the panel vibration Increases 
considerably and critical frequency behavior Is evident around 625 Hz. 

It should be noted that SEA often overpredicts in the critical frequency 
region. This simple example Illustrates some of the possibilities of SEA 
with the help of this computer program. 

A more extensive example Is that of the MEA done In the previous phase. 
The input for this six-element model (Fig. 8-1) Is extensive and much too 
elaborate to be explained In detail here, but the breakdown of the sub- 
elements used is given In Table 3-3. The Input listing Is given In 
Appendix III and the output is shown In Table 3-4. 

A mechanical vibration Input may also be applied If the known vibration 
level 1$ Included as an additional element. Load sheet 2 card column 5 
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STATISTICAL ENERGY ANA.YSIS OF COMPLEX STRUCTURES 


RECORD 

V UMBER DATA READ FROM UNIT 3 


1 NUMBER OF ELEMENTS = 2 

NUMBER 3 F ANALYSIS ? REQJE NC IES = 11 
FIRST ANALYSIS - REOUENCY = 2.50000E*02 

TYPE OF UNITS = 

TYPE OF OUTPUT = PSD 

2 ELEMENT NUMBER = 1 

NUMBER OF SUb-ELEMEV’S = 2 

TYPE OF EXCITATION = A 
TYPE OF MECHANIOAL INPUT = 

DAMPING = l.OOOQOE-Qi 
SLOPE = -1.0 3000E+00 

STARTING FREQUENCY = 2.50000E*02 

SOUND PRESSURE LEVELS 

3 l.4?OOOE*32 1.4r5Q0E«-02 1.4b000F»J2 1.4«OOOE>02 

1 • 4950 0E*G2 l.5C00QE*02 1.50000E*02 1.50000E*02 

4 l.S0000E>02 1* 53000E+02 1.50000E*02 0. 

0. 0. 0. o. 

J SUB-ELEMENT NUMBER = 1 

T»PE 0 F SUB-ELEMENT = P 
OENSITY = 2.G17Q0E-04 

MODULUS O r ELASTICITY = 1.00000E*07 

THICKNESS = 4.00 030E-02 

AREA = 1. 2 ?600E* 03 

POISSONS RATIO = 3.30OC0E-OI 

LENGTH = 0. 

PRESSURE = 0. 

6 STIFFNESS DEDUCTION REOUIREO = F 

RADIUS = 0. 

VOLUME = 0. 

SPEED OF SOUND IN ROOM MEDIUM = 1.34000E-HJ4 

ADDED MASS = 1.30S00E-O2 

7 SUB-ELEMENT NUMBER = 2 

TYFE OF SUB-ELEMENT = P 
DENSITY = 2.61700E-04 

MODULUS O r ELASTICITY = 1.00000E*07 

THICKNESS = G.30000E-02 

AREA £ 5 . 0 H j 0 0 C ♦ 0 2 

POISSONS RATIO = 3.30000E-01 

LENGTH = 0. 

PRESSURE = 0. 

R STIFFNESS REDJCTIOV REQUIRED = F 

RAOIUS = 0. 

VOLUME - 3. 

SPEED OF SOUND IN ROOM MEOIUM = 0. 

ADOED MASS = 0. 

(Continued) 

Figure 3-11. Two-Element SEA Program Output - External Acoustic Excitation 
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10 


ELEMENT NUMBER - 2 

NUMBER OF SUB-ELEMENTS = 1 

TYPE OF EXCITATION = 

TYPE OF MECHANICAL INPUT = 
DAMPING = 5.00000E-02 

SLOPE = 0, 

STARTING FREQUENCY = 0. 

SJB-ELEMfc^T NUMBER = 1 

TYPE OF SJ8-ELEHEN* = P 
DENSITY = 3.49Q00E-05 

MODULUS Or ELASTICITY = 
THICKNESS = 3*00 010E-01 

AREA = 5.1000QF»02 


1*0000 UE^P 7 



POISSONS RAT TO 
LENGTH = 0. 

PRESSURE = 0. 

“ 3. 30000E-0 1 

11 

STIFFNESS REOJCTIOM REQUIRED = 
RADIUS = 0. 

VOLUME = 0. 


SPEED OF SOUND 

TN ROOM MEDIUM = 


ADJEO MASS = 4 

•1SQ00E-02 

CENTER 

MODAL DENSITY - MODE 

FREQIH7) 

ELEMENT 1 

E.EMEN' r 2 

250.00 

5. 3H653E-J 2 

R.26429E-04 

312.50 

5. 36653E-0 2 

9.2b' 29E-04 

400.00 

5.39653E-32 

R.26429E-04 

300.00 

5 • 30653E-0 2 

R • 26429E-04 

625.00 

5, 3F653E-0 2 

H • 2 6429E-04 

767.50 

5.3E653E-02 

6.26429E-04 

1000.00 

5. 3R653E-82 

9.26429E-04 

1250. 00 

5. 3R653E-0? 

9.26429E-04 

1575.00 

5.3e>653E-02 

R.26429E-04 

2000. 00 

5.36653E-02 

A . 26429E-04 

2500.00 

5. 33653E- 9 2 

R. 26 429E-G4 


RECORD 

NUMBER DATA READ FROM UNIT 3 


12 FIRS 7 ELEMENT i l 
SECOND ELEMENT = 2 

TYPE OF JOINT = BJ 
NLMBER OF SIDES = 0 

JOINT LENGTH = 3*60 OOOE *0 1 

THICKNESS OF FIRST ELtMENT = 4.00000E-02 

THICKNESS OF SECONO ELEMENT = 3.00000E-01 

ACOUSTIC SPACE OENSITY = I). 

BEAM LENGTH = 0. 

INSERTION LOSS FACTOR = 2.0OG00C*0Q 


Figure 3-11 (Continued) 
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Table 3-1 


TWO-ELEMENT SEA VIBRATION PREDICTION 
EXTERNAL ACOUSTIC EXCITATION 


Center 

Frequency 


PSD Levels (G 2 /Hz) 


(Hz) 

Element 1 

Element 2 

250.0 

4.90675E+01 

3.38593E-01 

312.5 

3.59925E+01 

2. 37220E-01 

400.0 

2.5 3599 E+01 

1. 58301 E-01 

500.0 

2.11177E+01 

1.25095E-01 

625.0 

1.53027E+01 

8.85583E-02 

787.5 

1.1 7827E+01 

6.21 389 E-02 

1000.0 

7. 83661 E+01 

3.86783E-02 

1250.0 

5.44536E+00 

2.51832E-02 

1575.0 

3.81675E+00 

1.64492E-02 

2000.0 

2.72180E+00 

1.08695E-02 

2500.0 

2.05294E+00 

7.61286E-03 
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for that element would contain an M and the vibration spectrum descrip- 
tion would follow. The only restriction is that this element may not 
also have an acoustic input as it will be eliminated in the solution 
because the energy level is already known. 

The possible permutations and combinations of elements, sub-elements, and 
other factors which this computer program can create go far beyond the 
ability to document in this report. These few examples give some insight 
into the processes involved in the performance of SEA using this program. 
Although determination of some of the parameters such as damping and 
insertion loss is still a difficult and often obscure process which 
requires substantial future improvement, this computer program provides a 
significant step toward streamlining and simplifying SEA for the analyst. 
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Table 3-2 

TWO-ELEMENT SEA VIBRATION PREDICTION 
EXTERNAL AND INTERNAL ACOUSTIC EXCITATION 


PSD Levels (G 2 /Hz) 
Element 1 


Center 

Frequency 

(Hz) 

250.0 

312.5 

400.0 

500.0 

625.0 

787.5 
1000.0 

1250.0 

1575.0 

2000.0 
2500.0 


4.96984E+01 
3.67160E+01 
2.62012E+01 
2.26604E+01 
2.05716E+01 
1.40047E+01 
8. 75443E+00 
5. 75084E+00 
3.94377E+00 
2. 76144E+00 
2.06547E+00 


Element 2 

9.69548E-01 
1.04475E+00 
1.21999E+00 
2.29959E+00 
7.59740E+00 
3.98634E+00 
1.86390E+00 
7.03965E-01 
3. 33035E-01 
1.22184E-01 
4.69292E-02 
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Table 3-3 

SUB-ELEMENT BREAKDOWN FOR MEA TEST 


ELEMENT 1 


Sub-element 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


} 



Thermal Panels 
Orbiter Interface Panels 


Support Structure 


Radiator 


Interface Support Structure 
Signal Distributor Panel 
Signal Distribution Box 


Support Brace Assembly 
and Gusset 

Pressure Sensor and Voltage 
Regulator Panels 


ELEMENT 2 
Sub-element 
1 

2 • 

3 > 

4 - 

ELEMENT 3 
Sub-element 
1 
2 
3 

ELEMENT 4 
Sub-element 
1 

ELEMENT 5 
Sub-element 


ELEMENT 6 
Sub-element 

n 


CASE 

Power Distribution Panel 
Power Distribution Box 

Data Acquisition Cold 
Plate 


Solenoid Panel 


Battery Cold Plate 


Experiment Mounting 
Plates 
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Table 3-4 

SEA COMPUTER PROGRAM OUTPUT 
FOR MEA/ ACOUSTIC EXCITATION CASE 


CENTER 


PSD LEVELS 

(G**27HZ) 


FREQCHZ ) 

ELEMENT 1 

element 2 

ELCMENT 3 

ELEMENT 4 

31.50 

1. 3711 7E + 0 1 

6.40344E-02 

1. 23606E-32 

7 • 743 95E-0 2 

35.38 

1.39434E*0i 

5.33728E-02 

i. 16254E-02 

7.03336E-02 

50.40 

l»35352E*3l 

5*0 1959E-02 

1. 03534E-02 

6.006 38E-02 

63.00 

1. 378 34E*01 

4.58 005E-02 

9. J5592E-33 

5.434 42E- 02 

78.75 

1 • 40 49 5E + 0 1 

4.182 26 £-02 

9.20447E-03 

4.908 58E-02 

5 9.23 

l. 41 14 4E *0 1 

3. 7 48 61 £- 02 

8. 5364°£-03 

4.33913E-02 

126* 0 0 

1. 39 8 7 3E ♦ 0 1 

3.30 1 IQE-02 

7. 79522E-03 

3.75635E-02 

157. 5C 

1.13790E+31 

2*404 71E-02 

5.87560E-03 

2.6*5926-02 

198.45 

1.0251 7E«U1 

1.9 3208E-02 

4.89243E-03 

2 . 1 11 04E-02 

252.00 

7 . 250BQE *00 

1.22133 £-02 

3* 20463E-03 

1. 3030 IE-02 

315.00 

4. 73Q43E + 00 

8.57629E-03 

2.19825E-02 

9.29518E-03 

3° 3. 7 5 

3.10091E + 00 

6 .0 48 40 E- 03 

1.51521E-03 

6.654 39E-03 

504.00 

1. 95 71 7E ♦ 0 0 

4.13B89E-03 

1 . Cl 10 7E -02 

4.62347E-03 

630.00 

8.211b 3E-01 

1.36311E-03 

4* 46 15 €E -0 4 

2. 11322E-03 

787.50 

3.47857E-01 

3.51231E-04 

1. 98 79 OE - 94 

9. 74010E-04 

c 52. 25 

9.28257E-02 

2.44964E-04 

5.E9C56E-05 

2.33312E-04 

1260.00 

4.96 OOOT -02 

1.41506E-04 

3. 15427E-05 

1.65277E-04 

1575.00 

2. 77460E-02 

3. 513 14 £-05 

1. 85 67 SE- 05 

1.3Q240E-04 

1984.50 

1.24869E-02 

4.13Q79E-05 

8.61 054E-06 

4.i5i 41E-05 

2520.00 

5. 7485 5E-0 3 

2.05J37E-05 

4.28 49 IE - 06 

2.453 36E-05 

3150. OC 

2.82 04 3E - 33 

1.06428E-05 

2.21222E-06 

1 .30057E-05 

3937. 50 

1. 48151E-03 

6. 12303 £ - Ofa 

1. 22 29 7E-06 

7. 37309E-06 

5040.00 

8.43835E-04 

3.7 7799E-06 

7. 28 C86E-07 

4.56355E-06 

CENTER 


f»S3 LEVELS 

<6 • * 2/HZ > 


FR LQ< HZ » 

ELEMENT 5 

E-EME'i T & 



31.5 0 

3.C0958E-03 

3.5 39 44 E- 03 



38.38 

3. 5584 IE- 0 1 

8.381 72E-03 



5C.40 

3. 13632E- 03 

7.80207E-03 



63.00 

2.9 2 061E-33 

7.24503E-03 



78.75 

2.7171 OE- 0 3 

6.722 58 E-03 



85.23 

2.47o60E-03 

6.1 1710E-03 



126.00 

2 • 21 90 IE -03 

5.46337E- 03 



157.50 

1.63920E-03 

4*02753£-03 



196.45 

1 • 33442E-03 

3.27234E-03 



252.00 

8. 53722E-04 

2.089 91 l -03 



315.00 

5. 94798E- 04 

1.45769E-03 



39 3. 75 

4. 161S9E-04 

1.02111E-03 



504.00 

2. 82 109E-04 

6.93161E-04 



630.00 

l. 26179E-04 

3.10443E-04 



787.50 

5.69414E-05 

1.40254E-04 



992.25 

1.62123E-05 

4.0 0 061 £ -05 



1260.00 

9.2555/E-06 

2.28779E-05 



1575.00 

5.50 34 IE-06 

1.36259E-05 



1984.50 

2.63620E-06 

6. 533 79E-06 



2520.00 

1.29 3 2 4E-06 

3.21405£*0b 



3150.00 

6.72/12E-07 

1.&7509E-06 



3937.50 

3. 74282E-07 

9.33860E-07 



5040.00 

2.26442E-0 7 

5.67554E-07 
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Section 4 

SEA COMPUTER PROGRAM DESCRIPTION AND USAGE 

This section describes the computer code and the input requirements for 
problem solutions. 

The program was written in ASCII Fortran V for the MS'C UN I VAC 1108 with 
the EXEC8 operating system resident in the computer. The program is 
designed to be run in the demand/interactive mode using the system editor 
to assemble and sequence the load sheet information. 

The SEA computer code uses a main program whose function is essentially 
to call the various subroutines in the correct order. First called are 
two subroutines, EPROP and JINPUT, which read user-supplied data. Then 
for each analysis frequency, three other subroutines (JPROP, EXCITE, and 
ANSWER) are called to perform the various calculations. The subroutine 
ANSWER calls the UNIVAC library subroutine GASSEM from SYS $*MATHSTAT$. 
Finally, subroutine RITER is called to print the solution. The complete 
program listing is given in Appendix I and the program flow charts are 
given in Appendix II. 

4.1 SUBROUTINE DESCRIPTION 

Subroutine EPROP reads user-supplied data giving element and sub-element 
properties and calls one of five subroutines (BEAM, MEMBR , PI ATE, ROOM, 
or CYLIN) to calculate element modal densities. See Appendix I for the 
list of variables read by EPROP. EPROP first reads a record giving the 
number of elements, number of analysis frequencies, the first frequency, 
and the type of units (metric or English). EPROP checks that the number 
of elements and analysis frequencies Is within the bounds allowed by the 
program, since memory Is reserved with maximum values In mind. Then EPROP 
creates a table of analysis frequencies by multiplying each element of a 
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predefined 1/3 octave table by the first frequency. The predefined table 
consists of a series of frequencies, each of which is approximately 
one-third octave higher than the preceding one, and the first one of 
which is equal to 1. EPROP then reads element properties. If the type 
of excitation Indicated on this read operation Is acoustic, EPROP then 
reads a list of sound pressure levels for each analysis frequency for 
this element. Next, sub-element properties are read. Included in sub- 
element properties is the type of sub-element. Depending on whether this 
is a beam, membrane, plate, room or cylinder, EPROP calls the appropriate 
subroutine to calculate the modal density. 

Subroutines BEAM, MEMBR, PLATE, ROOM and CYUN calculate the modal 
density for a sub-element which is a beam, membrane, plate, room or 
cylinder, respectively. First, the part of the equation which is not 
frequency dependent is calculated. If stiffness reduction has been 
indicated, this partial value is multiplied by 1/ 2. Then for each 
analysis frequency, the rest of the modal density equation is computed 
and the result summed to element modal density. The sub-element mass is 
summed to the element mass. If this is the first sub-element for the 
given element, it is assumed to be the main sub-element and element 
properties other than mass are set equal to the properties of this sub- 
element. The following equations are used to calculate modal densities: 



, . A / Ft 2 K ' 1/z 

Plate: n(») = 

Vui 2 

Room: n(«) • 

Cylindrical shell: 
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Symbol FORTRAN Name 


Description 


n 

L 

u> 

E 

t 

o 

A 

S 

Y 

V 
c 


N 

Modal Density 

L 

Length 

OMEGA 

2ft Times the Frequency 

E 

Modulus of Elasticity 

T 

Thickness 

RHO 

Density 

A 

Area 

S 

Pressure 

GAMMA 

Poisson's Ratio 

V 

Volume 

C 

Speed of Sound 


Subroutine JINPUT reads joint properties supplied by the user. See 
Appendix I for a list of the variables read by JINPUT. JINPUT checks 
that element numbers are within range (i.e., less than or equal to the 
number of elements in the system) and that no pair of element numbers is 
input more than once. JINPUT keeps a running total of the number of 
pairs input and continues reading until the end of file is encountered. 
JINPUT checks to see that at least one pair of elements was read. 


Subroutine JPROP calculates coupling coefficients based on the data read 
by JINPUT, using one of the following equations, dependinq on the type of 
joint as read by JINPUT: 


Plate to plate at right angles: 



Beam to plate (cantilevered): 

2irfb 

♦i* “ NIC 

Plate to acoustic space: 


/4. 33ttc?* ' 

W2£j£ \ 

Vfw a V 2 - 

H p A / 
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Bolted or riveted joints: 

1) Calculate PHI as if for a plate-to-plate rigid joint. 

2) Reduce PHI by insertion loss factor. 


Symbol 

FORTRAN Name 

Description 

4 

PHI 

Coupling Coefficient 

L 

Jl 

Jor>t Length 

N 

MODES 

Number of Modes in Bandwidth 

ID 

OMEGA 

2t Times the Frequency 

t 

T 

Thickness 

E 

EE 

Modulus of Elasticity 

p 

DENSE 

Density 

r 

E GAMMA 

Poissr 'vs Ratio 

i 

TAU 

Thi:<ness Ratio 

n 

N 

*'odal Density 

f 

FREQ 

Frequency 

b 

BW 

Beam Width 

c 

EC 

Speed of Sound in Room Medium 

V 

VOL 

Volume 

a 

NS 

Number of Sides 

a 

SIGMA 

Radiation Efficiency 

P A 

ASD 

Acoustic Space Density 

Function SIGF returns a value for 

thd radiation efficiency of a panel. 

single argument, X, 

is passed to SIGF. The value of X is the ratio of 


the analysis frequency to the critical frequency. An internal table of 
16 values of the log of the radiation efficiency for 0 < X < 4 is main- 
tained. The first value, 0, is the value of X for which SIGF is a 


minimum. When X » 0, SIGF ■ -1.8. Each subsequent value of the internal 
table is the value of X for which SIGF increases by 0.2 over the previous 
value. The 13th value in the table is 1. Since this is the l?th value 
after the lit, SIGF - 12 * 0.2 + (-1.8) - 0.6 when X - 1. This is the 
maximum value of SIGF. SIGF * 0 for X - 4, the final value of the table. 
The value of SIGF is calculated by finding the least value of the table 
that is greater than X. This value and the precedinq one give two values 
of SIGF that differ by 0.?. Linear interpolation is then used to find 
the actual value of SIGF. 

Subroutine EXCITE calculates acoustic and mechanical energy inputs for 
the elements of the system. These values are initially set to 0. Then 
for each element, the acoustic or mechanical energy input is calculated 
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according to the following equations, depending on the type of excitation 
as read by EPROP: 


Acoustic: s . .H 

u) m 

<pS - 10 SPL/10 (8.41 x 10' ia ) 


Mechanical in g rms : E -^g 2 r(ns g 2 


Mechanical in 

PSD Levels: E 

- I" PSD g 2 f 
u 2 4.33 

Symbol 

FORTRAN Name 

Description 

S 

S 

Acoustic Energy Input 

c 

EC 

Speed of Sound In Room Medium 

A 

A 

Surface Area Exposed to Sound Field 

0 

SIGMA 

Radiation Efficiency 

N 

MODES 

Number of Modes In Bandwidth 

u 

OMEGA 

2t Times the Frequency 

m 

MSUB 

Mass 

SPL 

SPL 

Sound Pressure Level 

E 

E 

Element Energy Level 

9rm$ 

MECH 

Mechanical Input 

g 

G 

Gravitational Constant 

PSD 

MECH 

Mechanical Input 

f 

FREQ 

Frequency 

Subroutine ANSWER solves the SEA system of equations for element energy 

levels. First, element damping is 

determined. If the damping is con- 

stant, it Is equal 
equation Is used: 

to the value read by EPROP. Otherwise the following 

v ■ "f(y s 

Symbol 

FORTRAN Name 

Description 

n 2 

ETA? 

Element Damping 

nf 

ETA 

Constant Level of Damping 

f s 

SFREQ 

Start Frequency 

s 

SLOPE 

Slope 


Next, the elements of the alpha matrix of equation 1, Section 1, are 
calculated. 
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If there are any elements for which the energy levels (the E matrix) are 
already known, these are eliminated from the equation as shown by the 
following example: Suppose the system has six elements and the second 

and fifth have known energy levels. Then the equation becomes: 


®L1 °13 a 14 °16 


N 


S 1 “ ° 1 2 E 2 ' ^5] 

°31 933 934 936 


E 3 


s 3 " 932^ ~ °^5 E 5 1 

9*1 °H3 9*4 9*6 


E 4 

^ 28 

s 4 9 9*5 E 5 

961 963 964 956 


E 6 J 


s 6 - 95^2 “ 965 e 5. 


Since (for example) S x - + + 0^3 ♦ ♦ a^Eg + a^Eg, it 

be seen that this has the same solution as the original equation. ANSWER 
recalculates the values of the S matrix and calls SOLVE to eliminate the 
unnecessary , ows and columns from the matrices and find the solution. 

The solution is then used to calculate the average acceleration with the 
formula 


Symbol 

FORTRAN Name 

Description 

<a> 

ABAR 

Average Acceleration 

E 

E 

Element Energy Levels 

O) 

OMEGA 

2t Times the Frequency 

m 

MASS 

Mass 


Subroutine SOLVE solves the SEA system of equations. It copies the alpha 
and S matrices to new matrices, leaving out those rows and columns which 
were to be eliminated. It then calls MATHSTAT library subroutine GASSEM 
to solve the equation. SOLVE then puts the solution In the element 
energy array E, and subroutine RITER prints the solution. 

4.2 PROGRAM USAGE 

4.2.1 Deck Setup and Sequence 

At present, the source program resides on element C of file SI, so that 
It Is necessary to compile and collect It before execution. The program 


can 


\ 
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reads the SEA load sheet data on logical unit 3, which must be created by 
the user with the text editor or data processor, using the Q option so 
that file 3 is In ASCII code. The following sequence of control state- 
ments illustrates the creation of file 3 and the excecution of the 
program: 


0RUN ... 

0AS6.C 3. 

0ED.IQ 3. 

Statements creating file 3 
EXIT 

9ASG.A SI 
0FTN.N Sl.C.REL 
0MAP.IN SYM.ABS 
LIB SYS$*MATHSTAT$. 

LIB SYS**MSFC$. 

LIB SYS$*MSFC$. 

END 

3XQT ABS 
9FIN 

4.2.2 Input - Drum/Disk 

The only input for the program Is on logical unit 3, which contains the 
ucer'; data. This file consists of the following five types of records: 

1. Initial information used to process the other records. 

2. Element properties. 

3. Sound pressure levels or mechanical inputs. 

4. Sub-element properties. 

5. Joint properties. 

The arrangement of these records and their data elements Is shown on the 
following five loadsheet pages (Figures 4-1 through 4-5). 
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LOA0SHEET (1): HEADER CARD 

• One required per case 


f 
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Figure 4-1. Loadsheet 1 



LOADSHEET (2): SEA ELEMENT PROPERTIES CARD 




LOADSHEET (3): .XCITATION SPECTRUM CARD 

* Must immediately follow element properties if LS#2, CC5 f blank 



Figure 4-3. Loadsheet 3 



LOADSHEET (4): SUB-ELEMENT PROPERTIES, CARDS 1 AND 2 

• Must be consistent with Loadsheet 2, CC 3-4 (SE z 1) 

• Two cards (records) per sub-element 
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LOADSHEET (5): JOINT PROPERTIES 
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Figure 4-5. Loadsheet 5 




4.2.3 C'Jtput - Printout 

The output consists of a list of input parameters, calculated modal den- 
sities, and calculated PSD levels (or G ), if this was specified on 

input by the user) for each element for each analysis frequency, arranged 
in five columns below a heading. The first column is the frequencies, 
the other four columns list the PSD levels for elements 1 to 4. This is 
followed by five more columns below a new heading and containing the fre- 
quencies and PSD levels for elements 5 to 8. This is repeated until all 
the elements have been listed. 

4.2.4 Program Diagnostic Messages 

The program can produce the error messages listed below. Suggestions for 
their cause and correction are also given. Lower case x's indicate val- 
ues (usually, but not always, integers) that depend on the particular 
error. 


*** ERROR *■** THE INPUT FILE IS EMPTY. THE END OF FILE WAS 
ENCOUNTERED WHILE TRYING TO READ THE FIRST INPUT RECORD. THIS 
ERROR WAS DISCOVERED BY SUBROUTINE EPROP. 

*** ERROR *** WHILE TRYING TO READ ELEMENT PROPERTIES, THE END OF 
FILE WAS ENCOUNTERED BEFORE INPUT RECORD xxxx. THIS ERROR WAS 
DISCOVERED BY SUBROUTINE EPROP. 

*** ERROR *** WHILE TRYING TO READ SUB-ELEMENT PROPERTIES, THE END 
OF FILE WAS ENCOUNTERED BEFORE INPUT RECORD xxxx. THIS ERROR WAS 
DISCOVERED BY SUBROUTINE EPROP. 

*** ERROR *** WHILE TRYING TO READ SOUND PRESSURE LEVELS, THE END 
OF FILE WAS ENCOUNTERED BEFORE INPUT RECORD xxxx. THIS ERROR WAS 
DISCOVERED BY SUBROUTINE EPROP. 

*** ERROR *** WHILE TRYING TO READ MECHANICAL INPUTS, THE END OF 
FILE WAS ENCOUNTERED 8EF0RE INPUT RECORD xxxx. THIS ERROR WAS 
DISCOVERED BY SUBROUTINE EPROP. 
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*** ERROR *** THE END OF FILE WAS REACHED BEFORE ANY INFORMATION 
ABOUT JOINT PROPERTIES WAS READ. THIS ERROR WAS DISCOVERED BY 
SUBROUTINE JINPUT. 


These errors are the result of an incomplete data file. The data file 
cannot end before joint properties for at least one pair of elements are 
read. The message indicates what kind of data the program was looking 
for when the end of file was encountered. 

*** ERROR *** A FORTRAN ERROR OCCURRED WHILE TRYING TO READ 
ELEMENT PROPERTIES ON INPUT RECORD xxxx. THIS ERROR WAS DISCOVERED 
BY SUBROUTINE EPROP. 

*** ERROR *** A FORTRAN ERROR OCCURRED WHILE TRYING TO READ 
SUB-ELEMENT PROPERTIES ON INPUT RECORES xxxx AND xxxx. THIS ERROR 
WAS DISCOVERED BY SUBROUTINE EPROP. 

*** ERROR *** A FORTRAN ERROR OCCURRED WHILE TRYING TO READ THE 
FIRST INPUT RECORD. THIS ERROR WAS DISCOVERED BY SUBROUTINE EPROP. 

*** ERROR *** A FORTRAN ERROR OCCURRED WHILE TRYING TO READ SOUND 
PRESSURE LEVELS ON OR BEFORE INPUT RECORD xxxx. THIS ERROR WAS 
DISCOVERED BY SUBROUTINE EPROP. 

*** ERROR *** A FORTRAN ERROR OCCURRED WHILE TRYING TO READ 
MECHANICAL INPUTS ON OR 8EF0RE INPUT RECORD xxxx. THIS ERROR WAS 
DISCOVERED BY SUBROUTINE EPROP. 

*** ERROR *** WHILE ATTEMPTING TO READ JOINT PROPERTIES, A FORTRAN 
ERROR OCCURRED ON INPUT RECORD xxxx. THIS ERROR WAS DISCOVERED BY 
SUBROUTINE JINPUT. 

A FORTRAN error is usually the result of invalid characters appearing in 
a data field; for example, alphabetic characters appearing where the 
program expects to read an integer. This is likely to occur if data 
RECORDS are missing or out of order; for example, if an element proper- 
ties record says there are four sub-elements for that element, but records 
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for only three sub-elements are present. A FORTRAN error on the first 
record may indicate that the data file was not given the name 3. If the 
external file name is not 3, the internal file name should be made 3 by : 
USE statement. 


*** ERROR *** THE NUMBER OF ELEMENTS WAS GIVEN AS xx. IT MUST BE 
BETWEEN 2 ANO 20 INCLUSIVE. THIS ERROR WAS DISCOVERED BY SUBROUTINE 
EPROP. 


*** ERROR *** THE NUMBER OF ANALYSIS FREQUENCIES WAS GIVEN AS xx. 
IT MUST BE BETWEEN 1 AND 40 INCLUSIVE. THIS ERROR WAS DISCOVERED 
BY SUBROUTINE EPROP. 

*** ERROR *** THE TYPE OF MECHANICAL INPUT GIVEN FOR ELEMENT xx ON 
INPUT RECORD xxxx IS xxx. THIS IS NOT A VALID TYPE. THE TYPE MUST 
BE RMS OR PSD. THIS ERROR WAS DISCOVERED BY SUBROUTINE EPROP. 

*** ERROR *** TYPE * x ON INPUT RECORD xxxx IS AN INVALID TYPE. 
TYPE MUST BE B, M, P, C, OR R. THIS ERROR WAS DISCOVERED BY 
SUBROUTINE EPROP. 

*** ERROR *** THE TYPE OF JOINT GIVEN FOR ELEMENT PAIR xx AND xx 
ON INPUT RECORD xxxx IS xx. THIS IS NOT A VALID TYPE. THE TYPE 
MUST BE PP, BP, BJ, OR PA. THIS ERROR WAS DISCOVERED BY SUBROUTINE 
JPROP. 


*** ERROR *** THE DETERMINANT OF THE SEA EQUATION MATRIX IS 0. 
HENCE THERE IS NO SOLUTION. THIS ERROR WILL CAUSE THE PROGRAM TO 
ABORT. THIS ERROR WAS DISCOVERED BY SUBROUTINE SOLVE. 

These messages are self explanatory. If any of the errors listed thus 
far occur, the subroutine in which they occur will continue processing. 
Before control is transferred to another subroutine, however, the follow- 
ing message is printed and the program is aborted through a CALL FERR 
statement: 

BECAUSE OF THE ERRORS LISTED ABOVE, THE SEA PROGRAM WILL ABORT. 
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The following errors will not cause the program to abort, but may put the 
results in error: 

*** WARNING *** ON INPUT RECORD xxxx THE ELEMENT NUMBER WAS GIVEN 
AS xx, WHICH IS OUT OF ORDER. IT HAS BEEN CHANGED TO xx. THIS 
ERROR WAS DISCOVERED BY SUBROUTINE EPROP. 

Element properties are stored in arrays in the order in which they are 
read. If elements are referenced in the joint properties section in any 
other order, that is, when they are referenced as one element of a pair, 
the results are almost certainly erroneous. 

*** WARNING *** ON INPUT RECORD xxxx THE SUB-ELEMENT NUMBER WAS 
GIVEN AS xx, WHICH IS OUT OF ORDER. IT SHOULD BE xx. THIS ERROR 
WAS DISCOVERED BY SUBROUTINE EPROP. 

Sub-element properties are not stored in arrays so this is not likely to 
result in an error. If other errors occur, this message probably indi- 
cates that there are missing or extraneous records. 

*** WARNING *** ON INPUT RECORD xxxx, ONE OR BOTH MEMBERS OF THE 
ELEMENT PAIR xx AND xx WAS EITHER LESS THAN 1 OR GREATER THAN xx, 
THE TOTAL NUMBER OF ELEMENTS. THIS PAIR WILL BE IGNORED. THIS 
ERROR WAS DISCOVERED BY SUBROUTINE JINPUT. 

*** WARNING *** THE ELEMENT PAIR xx AND xx ON INPUT RECORD xxxx 
WAS PREVIOUSLY READ ON INPUT RECORD xxxx. THE FIRST VALUES WILL BE 
USED. THIS ERROR WAS DISCOVERED BY SUBROUTINE JINPUT. 

*** WARNING *** ON INPUT RECORD xxxx, BOTH ELEMENT NUMBERS WERE 
GIVFN AS xx. THEY MUST BE DIFFERENT. THIS RECORD WILL BE IGNORED. 
THIS ERROR WAS DISCOVERED BY SUBROUTINE JINPUT. 

These errors result in the indicated record being ignored. Whether or 
not the results are erroneous depends on whether the indicated record is 
necessary to the results. 

In addition, data which causes overflow, negative arguments to the square 
root functions, etc., will produce ASCII FORTRAN diagnostics. 
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Appendix I 
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Appendix III 


SEA PROGRAM INPUT LIST 
MATERIALS EXPERIMENT ASSEMBLY - EXAMPLE 2 



original page is 

OF POOR QUALITY 


**-«-*- mi- ■ i W ru Cwm-m 


«.■ 'c-i. W J4>d, . 


STATISTICAL ENCRQY ANALYSIS OF COMPLEX STRUCTURES 


ACCORD 

NUMBER DATA READ FROM UNIT 3 


1 


— 2 


3 

4 

5 
4 


f 


0 


* 


NUMBER OF CLEMEN'S s 4 
NUMBER OF ANALYSIS FREQUENCIES s 23 
FIRST ANALYSIS FREQUENCY * 3.15000£*01 

TYPE OF UNITS = 

TYPE OF OUTPUT = PSD 
ELEMENT NUMBER = I 
NUMBER OF SUB-ELEMENTS = 20 
TYPE OF EXCITATION = A 
TYPE OF MECHANICAL INPUT * 

0AHPIM4 x 1.0000QE-02 
SLOPE s Q. 

STARTINQ FREQUENCY s 0. 

SOUND PRESSURE LEVELS 
I.2Q000E*02 1.22000E*02 l.24COCE*02 

1.28000E+02 1* 30 00 Q£*02 1'32COC£*02 

1« 3450 OE * 02 1# 3500 0E *02 U35C0CE*02 

!• 35000£*02 1.33000E*02 i*31C0C£*O2 

1» 24000E*02 1*25000C*02 U23C0CC*O2 

1* 19000£*02 !• 17000E* 02 1.15CO€C*02 

SUB-ELEMENT NUMBER * I 
TYPE OF SUB-ELEMENT * P 
OENSITT = 2.41700E-04 

MODULUS OF ELASTICITY x 1.00000C«0? 
THICKNESS * 3WOCOOE-02 
AREA = 2* 3547 0£*04 

POISSONS RATIO x 3.30000C-01 
LCNOTH « 0« 

PRESSURE * 0* 

STIFFNESS REDUCTION REQUIRED = F 
RAOIUS * 0. 

VOLUME * •• 

SPEED OF SOUNO IN ROOM MEDIUM = 1«3400«E*04 

ADDED MASS = 4.49000E-01 

SUB-ELEMENT NLNBER * 2 

TYPE OF SUB-ELEMENT * P 

DENSITY s 2«41700£-04 

MODULUS OF ELASTICITY x 1«OOOOQC*07 

THICKNESS * 1*90 OOOC-Ol 

AREA * 2« 14500C* 02 

POISSONS RATIC * 3«30000EMi 

LCNQTH * 0* 

PRESSURE « I. 

STIFFNESS REDUCTION REQUIRED * F 
RADIUS a | t 
VOLUME x 0. 

SPCCO OF SOUNO IN ROOM MCOlUM * C. 

AOOCO MASS « I. 

SUB -CLEMENT NUMBER * 3 

TYPE OF SUB-CLCNCNT « P 

DENSITY * 2«4170 0E-04 

MOOULUS OF ELASTICITY s 1.00000C*0! 

THICKNESS x 9«50 I00E-42 

AREA « 3« 3120OC* 03 

POISSONS RATIC x 3.30000C-01 


24000E*02 
330U0E *02 
3500 QE *02 
27 00 OE *02 
21000002 


10 




€ 

* 


11 


€ 


12 


« 


c 


13 



« 


19 


€ 


I 


M 


* J \ • •'«« j<kW »S* « — • 

LENSTH s 0. 

PRESSURE x 0* 

STIFFNESS REDUCTION RCOU1REO 3 F 
RADIUS s 0. 

VOLUME s 0. 

SPEED OF SOUND IN ROOM MEDIUM s 0* 
ADDED MASS 3 0. 

SUB-ELEMENT NUMBER = 9 

TYPE OF SUB-ELEMENT * P 

DENSITY s 2* 61 700E-04 

MOOULUS OF ELASTICITY s 1.00000E«07 

THICKNESS 3 2.50000E-01 

AREA 3 6*900001* 02 

POISSONS RATIC 3 3. 30000E-01 

LENSTH s B. 

PRESSURE 3 o. 

STIFFNESS REDUCTION REQUIRED 3 f 
RAOIUS 3 a* 

VOLUME s o. 

SPEED OF 9 ' UNO IN ROOM MEDIUM * (. 

AODEO MAS- 3 0. 

SUB-ELEMENT NUMBER 3 5 

TYPE OF SUB-ELEMENT > p 

OENSITY « 2.61700E-04 

MODULUS OF ELASTICITY s 1.00000E»0T 

THICKNESS 3 9* 30 00 OE -02 

AREA 3 3.930001*03 

POISSONS RATIC s 3.30000E-01 

LENSTH s a. 

PRESSURE « 0* 

STIFFNESS REDUCTION REOUIREO 3 f 
RAOIUS * 0. 

VOLUME 3 0. 

SPEEO OF SOUNO IN ROOM MEOIUH s (. 
AOOEO NASS 3 0. 

SUB-ELEMENT NUMBER > 6 

TYPE OF SUB-ELEMENT 3 p 

DENSITY 3 2.617001-09 

MOOULUS OF ELASTICITY s i.OOOOOE«OT 

THICKNESS 3 2.50000E-01 

AREA s 1.22813E*03 

POISSONS RATIO 3 3.30000E-01 

LENSTH 3 a. 

PRESSURE s 0. 

STIFFNESS REDUCTION REBJIREO 3 f 
RAOIUS 3 a. 

VOLUME 3 a. 

SPEEO OF SOUNO IN ROOM MEDIUM 3 (, 

AOOEO MASS 3 a. 

SUB-tLCNCNT NUMBER 3 7 

TYPE OF SUB-ELCMENT 3 p 

OENSITY 3 2.617S0E-B9 

MODULUS OF ELASTICITY s 1.0000QE*07 

THICKNESS 3 >.asaaaE-«i 

AREA 3 1.747201*03 

POISSONS RATIO 3 3. 300801-81 

LENSTH 3 a. 

PRESSURE 3 a. 

STIFFNESS REDUCTION REOUIREO 3 f 
RAOIUS • a. 

VOLUME 3 a. 

SPEEO OF SOUND IN ROOM NCOlUM 3 c. 
AOOEO NASS * a. 

SUB-ELEMENT NUMBER 3 a 
TYPE OF SUB-ELEMENT s p 
DENSITY a P.AWAAT-BA 





< 

€ 

« 

( 



) 0 ) '< 



NOOULUS Of CiASTICITT x 1.000QOC«01 

THICKNESS - 4*44 OQQC *02 

UCA » 2.044 '7C* 03 

POISSONS OAT IC * 3.3000QC-01 

LENSTH < • • 

PRESSURE x 0* 

STIFFNESS REDUCTION REQUIRED = f 
MAO XUS x 0* 

VOLUME * 0* 

SPEED Of SOUNO IN NOON MEDIUM * C. 
ADDED NASS * 0. 

SUB-ELEMENT NUMBER * 9 

Tff E Of SUB-ELEMENT * P 

DENSITY x 2.4170QE-04 

NOOULUS Of CtAST ICXTY s 1,0Q000C*0I 

THICKNESS x 1.2500QC-01 

AREA x 7«47320£+02 

POISSONS RATIC x 3.300Q0E-Q1 

LENSTH x 0. 

PRESSURE x 0* 

STIFFNESS REDUCTION REQUIRED s F 
RADIOS x 0. 

VOLONE x 0. 

SPCCO Of SOUND IN NOON MEOlUM = €• 

ADDED NASS x 0. 

SUB-CLEMENT NUMBER » 10 

TYPE Of SUB-CLEMENT * P 

DENSITY x 2.4170OC-O4 

MODULUS OF EiASTICITT = l.Q0000E«Q1 

THICKNESS x S.42000C-01 

AREA x 4. OOOOOL* 01 

POISSONS RATIC x 3.30000001 

LCN4TH x 0. 

PRESSURE « 0. 

STIFFNESS REDUCTION REQUIRED x F 
RADIUS x 0 • 

VOLUME x o. 

SPEEO OF SOUND IN ROOM MEDIUM = C. 
AOOEO NASS x 0. 

SUB-ELEMENT NUMBER = 11 

TYPE OF SUB-ELEMENT x P 

OENSITT x 2*41 700E -04 

NOOULUS OF ELASTICITY x 1.00000C«0! 

THICKNESS s 2*50 OOOE -01 

AREA x 1«45920E+Q3 

POISSONS RATIC x 3.30000E-0I 

LEN6TH x 4. 

PRESSURE x 0. 

ST IPPNCSS REDUCTION REQUIRED * F 
RADIUS x 0* 

VOLUME x •« 

SPECJ OF SOUND IN ROOM MEDIUM * C. 
AOOEO MASS x 0 • 

SUB-ELEMENT NUMBER x 12 

TYPE OF SUB-ELEMENT * P 

DENSITY x 2*41 10 0E-04 

NOOULUS Of ELASTICITY « 1.00000C«C) 

THICKNESS « 3.7SC00C-01 

ARC* x S.030e0C*01 

POISSONS RATIC x 3.300Q0C-01 

LENSTH x 0« 

PRESSURE x I, 

ST I FF NESS REDUCTION RCQJIRCO x F 
RAOIUS x 0« 

VOLUME * 0* 

AMrrn nr KniNn in moon urn ? un * c. 


I 


i 



AODED MASS a O. 

SUB-ELEMENT NUMB E» = 13 

TYPE Of SUB-CLEMENT * P 

OENSITY = 2*61700E-04 

MODULUS Of ELASTICITY * UOOOOQOO? 

THICKNESS * 1*90 00 OE *01 

AREA = 3* 54 1G0E* 02 

POISSONS RATIO a 3.3OO0OC-01 

LENGTH a 0. 

PRESSURE = 0. 

STIFFNESS REOUCTION REQUIRED a F 
RADIUS = 0. 

VOLUME a 0. 

SPEED OF SOUND IN ROOM REOXUM a C. 
ACOEO NASS a 0. 

SUB-ELEMENT NUMBER a 14 

TYPE OF SUB-CLEMENT a p 

OENSITY * 2*A1700C-04 

MODULUS OF ELASTICITY a l.OOOOOUQJ 

THICKNESS a 1.50000E-01 

AREA * 7.3A2OOE«02 

POISSONS RATIC = 3.30000E-01 

LENGTH a 0. 

PRESSURE « 0 • 

STIFFNESS REOUCTION 4E0JIRED a f 
RADIUS a 0 • 

VOLUME a 0. 

SPEED OF SOUNO IN ROOM MEDIUM a C. 
AOOCO MASS a 0* 

SUB-ELEMENT NUMBER a 10 

TYPE CF SUB-ELEMENT a p 

DENSITY a 2.A1700E-04 

MOOULUS OF ELASTICITY * U0O0OO£«0i 

THICKNESS a 2*OOOOOC-OJ 

AREA a 1«90OOOE*O2 

POISSONS RATIO a 3.3O0OOC-O1 

LENGTH a 0* 

PRESSURE a o. 

STIFFNESS REOUCTION REQUIRED a f 
RADIUS a 0. 

VOLUME a 0. 

SPEED OF SOUNO IN ROOM MEDIUM a ( # 
AOOCO MASS a 0. 

SUB-CLEMENT NUMBER a 10 

TYPE OF SU8~£L£NENT * P 

OENSITY a 2.GI700C-04 

MOOULUS OF ELASTICITY a 1»0OOOOE«O9 

THICKNESS a 2.50 00 OC -01 

AREA a U0710QC + 02 

PCISSONS RATIC > 3.30000C-01 

LENGTH a 0. 

PRESSURE a 0. 

STIFFNESS REOUCTION REQUIRED a f 
RAOIUS a 
VOLURE a Q. 

SPEED OF SOUNO IN ROOM MEDIUM a C. 
AODEO MASS a 0, 

SUB-ELEMENT NUMBER a 17 

TYPE OF SUB-ELEMENT a P 

OENSITY a 2.A17OOC-04 

MOOULUS OF ELASTICITY a UO0QQOE«OJ 

THICKNESS a 1.9O0OOC-01 

AREA < 7*915000 01 

POISSONS RAUO a 3.3OOO0E-O1 

LENGTH < 0« 

^klr r x n. 


QfWGr PAGE IS 
OF POOH QUALITY 


c 

c 

« 

c 


"•^1 • 





« 

f 


39 ST'FFNESS REDUCTION REQUIRED = F 

RADIUS = 0. 

VOLUME = 0* 

SPEED OF SOUND IN ROOM MEDIUM = C. 
ADDED MASS = 0. 

40 SUB-ELEMENT NUMBER = 13 

TYPE OF SUB-ELEMENT a P 
DENSITY = 2.61700E-04 

MODULUS OF ELASTICITY a 1.00000E*0) 

THICKNESS a 2.50000E-01 

AREA a 1* ■’1800E-01 

POISSCNS RATIO = 3.30900E-01 

LENGTH = 0. 

PRESSURE = 0. 

41 STIFFNESS REDUCTION REOUIREO a F 

RAOIUS = 0. 

VOLUME = 0. 

SPEED OF SOUND IN ROOM MEDIUM = C. 
ADOEO NASS a 0. 

42 SUB-ELEMENT NUMBER a 19 

TYPE OF SUB-ELEMEF" = P 
DENSITY = 2.61700E-04 

MODULUS OF ELASTICITY a l.bdOOOC«OI 

THICKNESS a 2.00000E-01 

AREA a 1.55250E»01 

POISSONS RATIO a 3.3000UE-01 

LENfiTH = 0. 

PRESSURE = 0. 

43 STIFFNESS REDUCTION REQUIRED = F 
RADIUS a 0. 

VOLUME = 0. 

SPEED OF SOUND IN ROOM MEDIUM a t. 
ADDEO NASS a 0. 

44 SUB-ELEMENT NUMBER a 2D 

TYPE OF SUB-ELEMENT a P 
OENSITY = 2.61700E-04 

MODULUS OF ELASTICITY a 1.00000E«07 
THICKNESS a 1.90000E-01 
AREA a 1.14474E*0J 
POISSONS RATIO a 3.30000E-01 
LENGTH a o. 

PRESSURE s 0« 

45 STIFFNESS REDUCTION REQUIRED a T 

RAOIUS a 0. 

VOLUME a Q. 

SPEED OF SOUND IN ROOM MEDIUM a £• 
ADDED MASS a 0. 

44 ELEMENT NLMBER a 2 

NUMBER OF SUB-ELEMENTS a 4 
TYPE OF EXCITATION a 

TYPE OF MECHANICAL INPUT s 

DAMPING a 1.00000E-02 
SLOPE a -B.304B0E-01 
STARTING FREQUENCY a 2.50000E«02 
47 v SUB-ELEMENT NLMBER a 1 

TYPE OF SUB-ELEMENT a p 
DENSITY a 2.41700E-04 
MODULUS OF ELASTICITY a 1«OOOOOE«OI 
THICKNESS a 1.90000C-01 
AREA a 3*S41GQE+02 
POISSONS RATIO = 3.30000E-01 

LENGTH = 0. 

PRESSURE a 0. 

46 STIFFNESS REDUCTION REQUIREO a F 

RADIUS a 0. 

VOI IIMF a n. 




SPEED OF SOUND IN ROOM MEDIUM = C. 
RODEO NASS a 8.60000E-03 

49 SUB-ELEMENT NUMBER = 2 
TTPE OF SUB-ELEMENT * P 
DENSITY s 2*6170 OE -04 

MODULUS OF ELASTICITY = 1.00000E«0» 

THICKNESS = 1.50000E-01 

AREA s 7. 3600 OE* 02 

POISSONS RATIC = 3.30000E-01 

LENGTH = 0. 

PRESSURE s 0. 

50 STIFFNESS REDUCTION REQUIRED = F 

RADIUS = 0. 

VOLUME a 0. 

SPEED OF SOUNO IN ROOM MEDIUM * C. 
AOOEO MASS a 0* 

51 SUB-ELEMENT NLMBCR a 3 
TYPE OF SUB-ELEMENT a P 
OCNSITY a 2*61 70 OE -04 

MODULUS OF ELASTICITY = 1»00000E«07 

THICKNESS - 2*00 000E-01 

AREA a 1.98000E*02 
POISSONS RATIO = 3.30000E-01 

LENGTH * 0* 

PRESSURE = 0* 

32 STIFFNESS REDUCTION REQUIRED = F 

RACIUS a 0. 

VOLUME = 0* 

SPEED OF SOUNO IN ROOM MEDIUM = C. 
ADDED NASS a 0. 

S3 SUB-ELEMENT NUMBER = 4 

TYPE OF SUB-ELEMENT a P 
DENSITY a 2*61 700E-04 
MODULUS OF ELASTICITY = 1.00000E«0> 

THICKNESS = 2. SO 00 OE -0 1 

AREA * 1* 07100E-02 

POISSONS RATIC = 3.30000E-01 

LENGTH - 0. 

PRESSURE =0. 

34 STIFFNESS REDUCTION REQUIREO = F 

RADIUS = 0. 

VOLUME * 0. 

SPECO OF SOUNO IN ROOM MEDIUM = C. 
AOOEO NASS a 0. 

55 ELEMENT NUMBER * 3 

NUMBER OF SUB-ELEMENTS a 3 
TYPE OF EXCITATION = 

TYPE OF MECHANICAL INPUT * 

DAMPING a I.OOOOOE-02 

SLOPE a -8.30480E-Q1 

STARTING FREQUENCY » 2.500006*82 

56 SUB-ELEMENT NUMBER = 1 
TYPE OF SUB-ELEMENT a p 
DENSITY a 2.61 700E-04 

HOOULUS OF ELASTICITY * 1.00000E*0? 

THICKNESS * 6.30000E-01 

AREA = B * T6300E*01 

POISSONS RATIC a 3.30000E-0L 

LENGTH s 0. 

PRESSURE * 0* 

37 STIFFNESS REDUCTION REQUIRED = T 

RADIUS a 0. 

VOLUME a 0. 

SPCCO OF SOUNO IN ROOM MEDIUM * C. 
AOOEO HASS « 1.44800E-01 

4* SIlR-FI F MPNT NLMMTR * ? 






55 


60 


61 


62 


63 


64 


65 


66 


61 


TYPE OF SUB-ELEMENT = P 

DENSITY = 2* 617 00E-04 

f’OOULUS OF ELASTICITY = 1.0000QE + 0? 

THICKNESS = 3*50 CQ0E-01 

AREA = 9* 16 700E* 01 

POISSONS RATIC = 3»30000E-01 

LENGTH = 0. 

PRESSURE s 0. 

STIFFNESS REDUCTION REQJIRE3 = T 
RADIUS - 0* 

VOLUME = 0. 

SPEED OF SOUND IN ROOM MEDIUM - i. 
ADDEO MASS = 0* 

SUB-ELEMENT NUMBER = 3 

TYPE OF SUB-ELEMENT = P 

OENSITY = 2.61700E-04 

MODULUS OF ELASTICITY n 1*0Q000E«07 

THICKNESS = 7* 00 00 OE -01 

AREA = 6* 219 20£> 02 

POISSONS RATIO = 3.300HOE-01 

LENGTH - 0. 

PRESSURE s 0* 

STIFFNESS REDUCTION REQUIRED = T 
RAOIUS * 0. 

VOLUME - 0. 

SPEED OF SOUND IN ROOM MEDIUM - C. 
ADDED MASS = 0, 

ELEMENT NUMBER * 4 

NUMBER OF SUB-ELEMENTS = 1 

TYPE OF EXCITATION = 

TYPt OF MECHANICAL INPUT * 

DAMPIN6 - 1 •OOOOQE-02 

SLOPE = -8*3048 0E-01 
STARTING FREQUENCY = 2.50Q00E*02 

SUB-ELEMENT NUMBER = 1 

TYPE OF SUB-ELEMENT = P 
DENSITY = 2.6I700E-04 

MODULUS OF ELASTICITY £ 1*0OOOOE«O? 
THICKNESS- 2*50 000E-01 
AREA - 7*0 5 440E* 02 

POISSONS RATIC = 3.30000E-01 

LENGTH = 0« 

PRESSURE = 0* 

STIFFNESS REDUCTION REQUIRED - T 
RAOIUS = 0* 

VOLUME - 0« 

SPEED OF SOUND IN ROOM MEDIUM - C. 
ADDEO MASS = 2.27000E-02 

ELEMENT NUNBER * 5 

NUMBER OF SUB-ELEMENTS = 3 

TYPE OF EXCITATION - 
TYPE OF MECHANICAL INPUT = 

DAMPING = 1 • 000Q0E-02 

SLOPE s -8# 3048 QE-01 
STARTING FREQUENCY - 2.50000£*Q2 

SUB-ELEMENT NLMBER - 1 

TYPE OF SUB-ELEMENT s P 
OENSITY = 2*61 700E-04 

MODULUS OF ELASTICITY = 1.00000E*07 

THICKNESS = 6*50 OOOE -01 

AREA = 8 • f 6300E ♦ 01 

POISSONS RATIC = 3.30000E-01 

LENGTH S 0* 

PRESSURE = 0* 

STIFFNESS REDUCTION REQUIRED = T 
RAnTtic r n* 
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VOLUME = 0. 

SPEED OF SOUND IN BOON MEDIUM s (• 

BODED MASS = l.3S790E*00 

SUB-ELEMENT NUMBER s 2 

TYPE OF SUB-ELEMENT » P 

DENSITY s 2.61700E-04 

MOOULUS OF ELASTICITY * i«Q0Q0QE«09 

THICKNESS = 3*50 COOE-Ol 

AREA s 1 «09960E*02 

POISSONS RATIO = 3.3000QE-01 

LENGTH 3 0. 

PRESSURE « *« 

STIFFNESS REOUCTION REOUIREO = T 
RAOIUS = 0. 

VOLUME « 0* 

SPEED OF SOUND IN ROOM NCOlUM = C. 
ADDED MASS 3 0. 

SUB-ELEMENT NUMBER 3 3 

TYPE OF SUB-ELEMENT * P 

DENSITY 3 2*61 700E-0A 

MODULUS OF ELASTICITY * UOOOOOC+OY 

THICKNESS 3 I • 00 00 OE -01 

AREA 3 6 * Ml I0E*02 

POISSCNS RATIO 3 3.30000E-0! 

LENGTH = 0. 

PRESSURE 3 0 . 

STIFFNESS REOUCTION RCOJIRER 3 T 
RAOIUS 3 0. 

VOLUME 3 0. 

SPEED OF SOUND IN ROOM MEDIUM * C. 
ADDED MASS = 0. 

ELEMENT NUMBER * 6 

NLMHtR OF SUB-ELEMENTS * 2 

TYPE OF EXCITATION * 

TYPE OF MECHANICAL INPUT s 
DAMPING r I •OOOOOE-02 
SLOPE 3 -d * 30 AR 0E-01 
STARTING FREQUENCY 3 2*50 000E + 02 

SUB-ELEMENT NUMBER » l 
TYPE OF SUB-ELEMENT x P 
DENSITY s 2*6170 OE -OR 
MOOULUS OF ELASTICITY • l.OOOOOC'G? 
THICKNESS s 5*00 00 OE -01 
AREA = ?• 2bd 00E ♦ 02 

POISSONS RATIC » 3*30060£-01 

LENGTH x 0* 

PRESSURE £ 0* 

STIFFNESS REDUCTION REOUIREO * T 
RAOIUS £ 0. 

VOLUME 3 o, 

SPEED OF SOUND IN ROOM NEOIUM j> C* 

AOOEO MASS £ i*30320t*00 

SUB-ELEMENT NLM8E R 3 2 

TYPE OF SUB-ELEMENT « P 

OENSITY s 2a 61 TOOL -OR 

MODULUS OF ELASTICITY « 1*OOOOOE«QT 

THICKNESS * 3* 75 COOL -01 

AREA £ T.26400E+02 

POISSONS RATIO 3 3.30000E-01 

LENGTH s 0. 

PRESSURE s 0* 

STIFFNESS REOUCTION RCQUIREO * T 
RAOIUS 3 0* 

VOLUME * 0* 

SPEEO OF SOUND IN ROOM MEOIUM » C, 

Arnrr> mass 3 fl* 
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CENTCR 
FREOlH/l 
31 . 50 
39.30 

50.00 

63.00 
78.15 
99.23 

126.00 

157.50 
196.95 

252.00 
315.10 
393.75 

500.00 

630.00 

787.50 
992.25 

1260.00 

1575.00 
i«e*.5o 

2520.00 

3150.00 
3537.50 

5040.00 


NODAL DENSITY - MOOES/ CRAO/SEO 


ELEMENT 1 
1. OT89 2E *00 
1.0B*92C«00 
1.08892E*00 
1.0Bh9?E*00 
1.0089 2E*00 
1.08892E*00 
1.05e92E*00 
1.0*892E»00 
1«0889 2£«00 
1.08892E*00 
1.08S92E«00 
1.08892E*00 
1. 0889 2E ♦00 
1.06392E.00 
1. 0869 2E*0 0 
1.0889 2C*0C 
1. 0889 2E*00 
1.03892E*00 
1 • 08892E *00 
1.08892E*00 
1.08 89 2E*00 
1.0do92E*00 
1.08892E*00 


ELEMENT 2 
1.09014E-02 
1.09014E-02 
1.09014E-02 
1.09014E-02 
1.09014E-02 
1.09014E-02 
1.09014E-02 
1.09014E-02 
1.09 014E-02 
1.09014E-02 
1.09014E-02 
1.09014E-02 
1.09014E-02 
1.09014E-02 
1.09014E-02 
1.09014E-02 
1.09014E-02 
1.09014E-02 
1.09014E-02 
1.09014E-02 
1.09014E-02 
I.Q9014E-02 
1.09014E-02 


ELEMENT 3 
l,i097!E-0; 
1 . J OS 75E-03 
1.50S7EE-03 
1.J0S75E-03 
l.;097!E-OI 
1.J0S7TE-C3 
1.I0S75E-C3 
1. .0975E-Q3 
1. 10975E«OI 
1.S0S7EE-C3 
1.J097IE-03 
1. £0°75E-03 
1. JOS 7TE-03 
1. JOS 75E-03 
1. i0S7!E-03 
1.J0S75E-O3 
1. J0S7IE-03 
l.J0«7«E-03 
1 • I097SE -0 3 
1. S0975E-03 
1.J0S73E-03 

i.josiee-o; 

1.JGS7IE-03 


CENTER NODAL DENSITY - NODE S/ IR AD/SEC ) 


FRtOCHZ) 

ELEMENT 5 

ELEMENT 8 

31.50 

1.32509E-03 

3.19268E-03 

39.38 

1.32509E-03 

3.19286E-03 

50. 0C 

1 • 32509E-0 3 

3.19288E-03 

63.00 

1.32509t-03 

3.19286E-03 

78.75 

1.32509E-03 

3.19288E>03 

99.23 

1. 32509E-03 

3.19286E-03 

126.00 

1. 32509E-03 

3.19288E-03 

157.5© 

1. 32509E-03 

3.19288E-03 

196.05 

1. 32509E-03 

3.19288E-03 

252.00 

1.32509E-03 

3.19286E-03 

315.CC 

1.32509E-03 

3.19288E-03 

393.75 

1.32509E-03 

3.19286E-0:i 

500.00 

1.32509E-03 

3.19288E-03 

63C.OO 

1. 32509E-03 

3.19286E-03 

787.50 

1.32509E-03 

3.19288E-03 

952.25 

1.32509E-03 

3.19268E-03 

1260. OC 

1.32509E-03 

3.19268E-03 

1575.00 

1.32509E-03 

3.19286E-03 

15P0.5C 

1.32509E-03 

3.15288E-03 

2520.00 

1.32509E-03 

3.19286E-03 

3150.00 

1.32509E-03 

3.19288E-03 

3937. 5C 

1.32509E-03 

3.19288E-03 

5000.00 

1.32509E-03 

3.19288C-03 


« 

f 

< 


RECORD 

NUMBER 


DATA READ FR3N UNIT J 


IT FIRST ELEMENT s 1 
SECOND ELEMENT = 2 

TYPE OF JOINT * BJ 
mi mbf» of stnrs = o 


ELEMCN T 4 
2.65614E-03 
2.65&14E-01 
2.65614E-03 
2.65614E-03 
2.65614E-0J 
2.65614E-03 
2.65'il4E-03 
2.65614E-03 
2. 656 14E-03 
2.65ol4E-03 
2.85614E-03 
2. 65614E-03 
2.65614E-03 
2.85614E-03 
2.65614E-03 
2. 65614E-03 
2.65814E-03 
2.85&14E-03 
2.65614E-03 
2.85814E-03 
2.65614E-0J 
2.65614E-03 
2.45814E-03 




JOINT LENGTH x 7.06TOOE*01 
THICKNESS OF FIRST ELEHENT = 2.5OC00E-01 

THICKNESS OF SECOND ELEMENT x 1.90000E-C1 
ACOUSTIC SPACE OENSITT = 0. 

BEAM LENGTH = 0. 

INSERTION LOSS FACTOR = 1.00000E*02 

76 FIRST CLEMENT x I 
SECOND ELEMENT x 3 
TYPE OF JOINT x BJ 
NUMBER OF SIDES = 0 

JOINT LENGTH = 1.23S40E*02 

THICKNESS OF FIRST ELEMENT s 2.50000E-0! 
THICKNESS OF SECOND ELEMENT x 7.0000CI-01 
ACOUSTIC SPACE DENSITY x o. 

SCAM LENGTH s 0. 

INSERTION LOSS FACTOR x 1.00000E*02 
n FIRST ELEMENT x 1 

SECOND ELEMENT s 4 
TYPE OF JOINT s BJ 
NLMBER OF SIDES = 0 

JOINT LENGTH x 4.00000E«00 
THICKNESS OF FIRST ELEMENT s 2»50000E*0! 
THICKNESS OF SECOND CLEMENT : 2.50000C-CI 

ACOUSTIC SPACE OENSITT s 0. 

BEAM LENGTH = 0. 

INSERTION LOSS FACTOR x 1.00000E*02 
60 FIRST ELEMENT x 4 

SECOND CLEMENT = 3 

TYPE OF JOINT x PP 
NUMBER OF SIDES = 0 

JOINT LENGTH = 1.6?600E*01 

THICKNESS OF FIRST ELEMENT = 7.00000C *0 1 

THICKNESS OF SECOND ELEMENT x 2.500001-01 
ACOUSTIC SPACE DENSITY x 0» 

BEAM LENGTH = 0. 

INSERTION LOSS FACTOR s 0. 

01 FIRST CLEMENT x 1 

SECOND CLEMENT x 5 
TYPE OF JOINT x BJ 
NLMRER OF SIDES x 0 
JOINT LENGTH x 1.08840E»02 
THICKNESS OF FIRST CLEMENT « 2#50000E-01 

thickness of second element s t.qooooe-oi 

ACOUSTIC SPACE OENSITT x 0. 

BEAM LENGTH x 0. 

INSERTION LOSS FACTOR = 1.00000E*02 

62 FIRST CLEMENT x 1 

SECONO ELEMENT » G 
TYPE OF JOINT x BJ 
NUMBEM OF SIDES = 0 

JOINT LENGTH s 2.G7500E*02 
THICKNESS OF FIHST ELEMENT s 2.50000E-01 
THICKNESS OF SECONO CLEMENT x 3.450001*01 
ACOUSTIC SPACE OENSITT x 0* 

BEAM LENGTH x 0. 

INSERTION LOSS FACTOR x 1.00000C*02 


« 



